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Abstract. In this article, we derive a rigorous characterization of the 
boundary-layer and interior separations in the Taylor-Couette-Poiseuille 
(TCP) flow. The results obtained provide a rigorous characterization 
on how, when and where the propagating Taylor vortices (PTV) are 
generated. In particular, contrary to what is commonly believed, we 
show that the PTV do not appear after the first dynamical bifurcation, 
and they appear only when the Taylor number is further increased to 
cross another critical value so that a structural bifurcation occurs. This 
structural bifurcation corresponds to the boundary-layer and interior 
separations of the flow structure in the physical space. 



1. Introduction 

Consider a viscous fluid between two coaxial rotating cylinders. The base 
(Couette) flow becomes unstable as soon as the rotation speed of the inner 
cylinder exceeds a critical value. This instability gives rise to a stationary 
axisymmetric counter-rotating vortices that fill the whole annular region. 
The associated flow is referred to as Taylor-Couette (TC) flow. When a 
through flow driven by a pressure gradient along the rotating axis is added, 
the resulting system can exhibit both convective and absolute instabilities. 
The base flow consists of a superposition of circular Couette flow and annular 
Poiseuille flow, called Couette-Poiseuille (CP) flow. The axial through- 
flow suppresses the basic stationary instability, and as the axial pressure 
gradient increases while the rotation speed of the inner cylinder is held 
fixed, the first bifurcation gives rise to a traveling train of axisymmetric 
Taylor vortices, commonly referred to as propagating Taylor vortices (PTV). 
Henceforth, the term Taylor-Couette-Poiseuille (TCP) flow is used to refer 
to all hydrodynamic phenomena pertaining to the open system described 
above; see among others Raguin and Georgiadis [13] and the references 
therein. 

In this article we rigorously characterize the stability and transitions of 
the CP and TCP flows in both the physical and phase spaces. The main 
focus is on a rigorous characterization on how, when and where the PTV 
are generated. 
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We first examine the existence, explicit formula, and basin of attractions 
of the secondary flows. The existence of a bifurcation to a steady state so- 
lution is classical, and can be proved by the classical Krasnoselskii theorem. 
The new ingredients here are the4 explicit formula and the basin of attrac- 
tion of the bifurcated solutions. These enable us to determine, in the next 
two sections, the asymptotic structure in the physical space of the solutions 
of the problem, leading to justification of the flow structure of the TCP flow. 

The main result of this article is on the rigorous characterization of the 
boundary-layer and interior separations, associated with the TCP flows. In 
particular, contrary to what is commonly believed, we show that the PTV 
do not appear after the first dynamical bifurcation, and they appear only 
when the Taylor number is further increased to cross another critical value so 
that a structural bifurcation occurs. This structural bifurcation corresponds 
to the boundary layer and interior separations of the flow structure in the 
physical space. This study gives another example linking the dynamics of 
fluid flows to the structure and its transitions in the physical spaces. 

The analysis is based on a geometric theory of two-dimensional (2D) in- 
compressible flows, and a bifurcation theory for nonlinear partial differential 
equations and both developed recently by the authors; see respectively [TT] 
and [12] and the references therein. 

The geometric theory of 2D incompressible flows was initiated by the au- 
thors to study the structure and its stability and transitions of 2-D incom- 
pressible fluid flows in the physical spaces. This program of study consists 
of research in two directions: 1) the study of the structure and its tran- 
sitions/evolutions of divergence- free vector fields, and 2) the study of the 
structure and its transitions of the flow fields of fluid flows governed by the 
Navier-Stokes equations or the Euler equations. The study in the first di- 
rection is more kinematic in nature, and the results and methods developed 
can naturally be applied to other problems of mathematical physics involv- 
ing divergence- free vector fields. In fluid dynamics context, the study in 
the second direction involves specific connections between the solutions of 
the Navier-Stokes or the Euler equations and flow structure in the physi- 
cal space. In other words, this area of research links the kinematics to the 
dynamics of fluid flows. This is unquestionably an important and difficult 
problem. Progresses have been made in several directions. First, a new 
rigorous characterization of boundary-layer separations for 2D viscous in- 
compressible flows is developed recently by the authors, in collaboration in 
part with Michael Ghil; see [12] and the references therein. Another example 
in this area is the structure (e.g. rolls) in the physical space in the Rayleigh- 
Benard convection, using the structural stability theorem developed in Area 
1) together with application of the aforementioned bifurcation theory; see 
[Hi [12] . We would like to mention that this article gives another example in 
this program of research. 

The dynamic bifurcation theory is centered at a new notion of bifurcation, 
called attractor bifurcation for dynamical systems, both finite dimensional 
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and infinite dimensional, together with new new strategies for the Lyapunov- 
Schmidt reduction and the center manifold reduction procedures. The bi- 
furcation theory has been applied to various problems from science and en- 
gineering, including, in particular, the Kuramoto-Sivashinshy equation, the 
Cahn-Hillard equation, the Ginzburg-Landau equation, Reaction-Diffusion 
equations in Biology and Chemistry, the Benard convection problem and 
the Taylor problem in classical fluid dynamics, and the some geophysical 
fluid dynamics problems. We mention the interested readers to a recent 
monograph by the authors [llj and the references therein. 

The paper is organized as follows. Section 2 introduces the physical prob- 
lem, TCP problem. Section 3 proves an abstract dynamical bifurcation 
theorem, which will be used in Section 4 to derive an explicit form and the 
basin of attractions of the bifurcated solutions. Sections 5 and 6 characterize 
the boundary-layer and interior separations associated with the TCP flow. 



2. COUETTE-POISEUILLE FLOW 



2.1. Governing equations. We consider the viscous flow between two ro- 
tating coaxial cylinders with an axial pressure gradient. The instability 
of this flow was first discussed by S. Goldstein in 1933 [5], and a further 
investigation was made by S. Chandrasekhar pQ. 

Let n and r 2 (n < r 2 ) be the radii of two coaxial cylinders, fii and 2 
the angular velocities with which the inner and the outer cylinders rotate 
respectively, and 



(2.1) 



p = VL 2 j^l\, T) = n/r 2 . 



The hydro dynamical equations governing an incompressible viscous fluid 
between two coaxial cylinders are the Navier-Stokes equations in the cylin- 
der polar coordinates (z,r,9), given by 



(2.2) 
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where v is the kinematic viscosity, p the density, u - 
field, p the pressure function, and 
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(u r ,ug,u z ) the velocity 
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In this article, we consider the case where a constant pressure gradient 
g| = po is applied in the z-direction, the cylinders rotate, and there are no 
radial motions. Under these conditions, (|2.2p admits a steady state solution: 



(2.3) U cp = (u z ,u r ,u e ) = (W,0,V), p = p(r,z), 
where 

(2.4) ifWVr, 

(2 5) i/ d f ^ + ^ V - 

dr \dr r J 

(2.6) „(' + iU W _B>, 

\dr r J dr p 



supplemented with the following boundary conditions 
(2.7) J (M ^ ) = (°'~^ 




0, W ) at r = n,r 2 , 



r 



at r = ri (i=l, 2). 

Here Wo > is a constant. 

We derive then from ([273]) — ([277]) that 

(2.8) U(r) = ( r?(l - M )- - (rf ~ H> 

(2.10) p(r, 2f) = p z + p y - F 2 (r)iir, 

where d (= r 2 — r\) is the gap width, p, and rj are as in (|2.ip . 
Hereafter, we always assume that 

(2.11) r? 2 > p, 

which is a necessary condition for the instability to occur. 

Classically, the flow described by U c = (0,0, V(r)) is called the Couette 
flow, and the flow by W p = (W(r),Q, 0) is called the Poiseuille flow. There- 
fore the flow described by (]2,3p is a super-position of the Couette flow and 
the Poiseuille flow, and is usually called the Couette-Poiseuille (CP) flow. 

The main objective of this article is to study the stability and transitions 
of the CP flow in both the physical and the phase spaces. 

For this purpose, we consider the perturbed state 

(2.12) W + u z , u r , V + u e , p + poz + p J ^V 2 dr. 
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Assume that the perturbation is axisymmetric and independent of 6, we 
derive from (|2.2ft that 



(2.13) 



9u z du z dW 1 dp 

— — + (u ■ \/)u z = vAu z - W — u r — — , 
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dz dr 
where u = (u z , u r ), and 



A - — -— — 

dr 2 r dr dz 2 ' 

,~ ^ 9 9 

The spatial domain for (|2.13p is M = (0,L) x (n,^), where L is the 
height of the cylinders. The physically sound boundary conditions are as 
follows 



(2.14) 



u = (u z ,u r ,u e ) = 
du r duo 



u z = 0, 



dz dz 



at r = t\ , T2 , 
at z = 0, L. 



2.2. Narrow-gap approximation. Consider the narrow-gap approxima- 
tion with p, > 0. Namely, we assume that the gap d = r 2 — r\ is small in 
comparison to the mean radius: 

(2.15) d = r 2 -r x < -(r 2 + ri). 

Under this assumption, we can neglect the terms having the factors r~ n 
(n > 1) in the equations. Let 



(2.16) 



1 — rf 
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By (l2~TTjh a > 0. We derive then from and flZ3]) that 



(2.17) V r/ + -y = -2a^i, 
(2,8) 

(2.19) W~-^_[(r-n)(r 2 -r) + Wb]. 



Replacing ug by y/aug in (|2.13p . we obtain the following approximation 
equations describing the flow between two cylinders with a narrow gap: 
(2.20) 

( du 

-— ^ + (u ■ V)u z =vAu z + - — (d - 2r)u r 
at 4pu 

Po (xxr . , -,sdu z 1 dp 

+ ~ A — {Wo + r(d-r))— — , 

Apv Oz p Oz 

du l (1 — LltT \ 

+ (u • V)u r =uAu r + 2y/a Q 1 I 1 - - — J Ug 

+ — {Wo + r(d-r))— — , 

Apv oz p Or 

+ (u ■ V)u e =uAu e + 2^fi lMf + p- (Wo + r(d - r)) ^ , 

divu = 0, 



where r = r — n, A = d 2 /dz 2 + d 2 /dr 2 . 

We need to consider the nondimensionalized form of (|2.20p . To this end, 

let 



x = x d (x = {z, r, r0)), 

t = t'd 2 /v, 

U = Uu/d (u = (u z ,U r ,Ug)), 

p = p'pu 2 /d 2 , 

Wo = Wod 2 . 
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Omitting the primes, the equations (I2.20p can be rewritten as 
( du. 



(2.21) 



at 

du r 
~dt 



■ Au z + 7(1 — 2r)u r 
+ 7(^0 + r(l - r))^j - g - (5 • V)u z , 



=Au r + A(l — (1 — fJ>)r)ug 

- 7 (W + r(l - r))^ - ^ - (u ■ V)ur, 
oz Or 

^ =Au e + \u r + 7 (Ty + r(l - r))^ - (u ■ V)u e , 
ot oz 

div u = 0, 



where A = vT. Here the Taylor number T and the nondimensional param- 
eter 7 are given by 



(2.22) 
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We remark here that the nondimensional parameter 7, proportional to the 
Reynolds number Re squared times the small gap d, is a small parameter. 

When the gap d = ri — t\ is small in comparison with the mean radius, we 
need not to distinguish the two equations (|2.13p and (|2.2ip to discuss their 
dynamic bifurcation. Therefore, we always consider the problem (|2.2ip with 
(gJID instead of (pJ3|) with (|234~j) . The equations ([231]) are supplemented 
with the following initial condition 



(2.23) 



u = ip at t = 0. 



3. Dynamic Bifurcation for Perturbed Systems 

3.1. Preliminaries. Let H\ and H be two Hilbert spaces, and H% — * H 
be a dense and compact inclusion. We consider the following nonlinear 
evolution equation: 




— = Au + XBu + C\u + G(u, A), 
u(0) = if, 



where u : [0, 00) — > H is the unknown function, A G R and e € M m (m > 
1) are the system parameters, and A, B : H\ — > .ff are linear operators, 
£| : i^i — > are parameterized linear operators continuously depending on 
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A E R, £ E R m , which satisfy 

' A : H\ — > H a sectorial operator, 
(3.2) < A a : H a — ► H the fractional power operators for a E 
H a = D(A a ) the domain of A a with H = H; 
f C\, B : Hq — > il bounded for some < 1, 



(3.3) 



£y^0 if e->0, VAe 



We know that H a C .ffg are compact inclusions for all a > f3, and the 
operators A + XB + £^ : i2i — ► are sectorial operators. Furthermore, 
we assume that the nonlinear terms G(-,A) : H a — > H for some a < 1 
are a family of parameterized C r bounded operators (r > 1), depending 
continuously on the parameter A E R, such that 

(3.4) G(u,X)=o(\\u\\ Ha ), VAeR. 

Let {S^li^o be an operator semigroup generated by (|3.ip . which enjoys 
the semigroup properties. Then the solution of (|3.1|) can be expressed as 
u(t, if) = S(t)ip for any t > 0. 

Definition 3.1. A sei T, C H is called an invariant set of Ii3.1\) ifS(t)Ti = £ 
/or any t > 0. ^4n invariant set T, C H of \3. 1\) is called an attractor if £ 
is compact, and there exists a neighborhood U C H of ' £ suc/i i/iaf for any 
ip £ U we have 

lim dist(S(t)(p, £) = in i? norm, 

t^oo 

and we say that £ attracts U . Furthermore, if U = H , then £ is called a 
global attractor of 113.1]) . 



A number /3 = /?i + i/3 E C is called an eigenvalue of a linear operator 
L : .f/i — ► H if there exist x,y £ Hi with i/O such that 

Lz = {z = x + iy). 

The space 

#0= (J {x,yE#i | (L-/3) n z = 0, z = x + iy} 

nGN 

is called the eigenspace of L, and x,y £ Ep are called eigenvectors of L. 

Definition 3.2. A linear mapping L* : Hi H is called the conjugate 
operator of L : Hi — ► , i/ 

(Lx,y) H = (x,L*y) H , V x,y e H x . 
A linear operator L : #1 H is called symmetric if L = L* . 
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3.2. Eigenvalue problem. Here we consider the eigenvalue problem for 
the perturbed linear operators 



Lx + C% 



L X = A + XB. 

Let the eigenvalues (counting multiplicities) of L\ and L\ be given respec- 
tively by 

{/3 fc (A) | A; = 1,2,...} CC, 

{^(A)|fe = l I 2,...}cC. 

Suppose that the first eigenvalue /?i(A) of L\ is simple near A = Ao, and 
the eigenvalues /?fc(A) (k = 1,2,...) satisfy the principle of exchange of 
stabilities (PES): 

< 
= 



(3.5) 



A (A) 



if A < A , 
if A = A , 



> if A > A 



o- 



Re(3 j (X )<0 Vj>2. 



(3.6) 

Let e\,e* x £ Hi be the eigenvectors of L\ and L* x corresponding to an eigen- 
value /3(A) respectively, i.e., 

[ Ae x + \Be x 



(3.7) 

We also assume that 
(3.8) 



A*e* x + XB*e* x 



/?(A)e A , 
= 0(A)e£. 

= 1. 



The following theorem is important for our discussion later, which pro- 
vides a differential formula for simple eigenvalues of linear operators. 



Theorem 3.3. Let the eigenvalue /3(A) in (3.1) be simple. Then /3(A) is 
differentiable at X, and with the assumption \3. 8\) we have 

(3.9) P'(X) = (Be x ,e* x ) H . 

Proof. By the genericity of simple eigenvalues of linear operators (see Kato 
[B] and Ma and Wang [H]), there is a number 5o > such that as \S\ < So, 
/3(A + 5) is also a simple eigenvalue of L\ + $. Let 

' Au + (X + S)Bu = 0(X + 5)u, 

u = e x + v s , 

J\vs\\h — >0 if S^O. 

We infer from <^ and (l3TT0|) that 

({A + XB)e\, e\) H + ((A + XB)v 5 , e\) H + S(Bu, e\) H 

= (i{X) + f3{X){vs,e\) H + S{Bu,e\) H 

= P(X + S)+f3(X + S)(v 5 ,e* x ) H . 



(3.10) 
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Thus the last equality implies that 
/3(\ + 8)-P{\) 



1 + (vs, e* x ) H ) = (B(e\ + v s ), e* x ) H . 



Then the theorem follows from (|3. 101) . 



□ 



Remark 3.4. For general linear operators L x : Hi — * H, if L x are differ- 



entiable on A with C x = ^ L x : Hi —> H bounded, then Theorem 
true as well, and 

(3.11) 



P'(X) = (C x e x ,e* x ) H . 



holds 



□ 



By Theorem 13.31 or (13. lip , for the perturbed linear operators L £ x we im- 
mediately obtain the following corollary. 

Corollary 3.5. Assume the PES (EO)) and ( TO)) , and ^ C\ exists. Then 
for each \e\ < 5 for some 5 > 0, there exists Aq, with lim^o A§ = Ao, such 
that the eigenvalues /3|(A) of L e x at Ag satisfy the following PES: 

'<0 if\<\ £ , 

(3.12) Pt(X)l =0 if\ = \ £ , 

>0 if\>\ £ , 

(3.13) Re^(\ E ) < 0, V 2 < j. 

Moreover, (5 £ (A) has the following expansion at A = Aq/ 
r/?f(A) = a e (A-A )+o(|A-A £ |), 
\a £ = ((B + ±£lc)e £ ,e* £ ) H , 

where e £ and e* satisfy 

(3.15) L%e e = 0, L £ x \e* s = 0. 



(3.14) 



3.3. Bifurcation for perturbed equations. We study now the dynamic 
bifurcation for the nonlinear evolution equation (|3.ip . which is a perturba- 
tion of the following equation 

du 

(3.16) — = Au + \Bu + G(u,\). 

dt 

Under the conditions (|3.7p and (|3.8p . let eo € -ffi be an eigenvector of L x 
at A = Aq: 



(3.17) 

We assume that 
(3.18) 



Acq + Ao-Beo = 0, 

f L x = A + \B 

G{u,X) 

(G(u,v),v)h = 



l e o| 



1. 



are symmetric, 
are bilinear, 
V u, v £ -Hq. 
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Here if G is bilinear, one can write G as G(-,-), which is linear for each 
argument. 

By Theorem 6.15 and Remark 6.8 in [11], together with Corollary 13.51 
above, we obtain immediately the following theorem. 

Theorem 3.6. Assume that \3.2\) ~ [3^\) and \3.1S\I hold true, and u = 
is a globally asymptotically stable equilibrium point of t3.16\) at X = Xq. 
Then there are constants 5i > and 62 > such that if |e| < 5\ and 
< A — Ag < 62, where Aq is given in A3.12\) and h3.13\) , then the following 
assertions hold true. 

(1) There exists an attractor = {u^, u^} C H of 113. where u\ and 
U2 are steady states given by 



(3.19) 



' u\ = fJi(A,e)e + vi(X,e), 
< U2 = -cr 2 (A, e)e + v 2 (X, e), 
k «i(A,e) =o(\o-i(X,£)\) €fli, i = l,2. 

Here crj(A,e) G M. can be expressed as 



y^( £ )+4/?f(A)T|fr(g)l , ,,,, |flen 
(3.20) { ai < 2 = 2C (\ b \>\Pi\)' 

b(e) = (G(e E ,X £ ),e* £ ) H , 

where e e ,e* are as in \3.15\) . /3f (A) as in 13.14% and C > is a 
constant. 

(2) attracts every bounded open set in H \ T, where T is the sta- 
ble manifold of u = with co-dimension one. Namely, H can be 
decomposed into two open sets U± and U^: 

h = u\ + U2, u± n U2 = 0, e r = du? n du£, 

with £ (i = 1,2), suc/i i/iai 

lim ip) - u^\\ H = 0, iftptUl {i= 1, 2), 

t— >oo 



where u(t,ip) is the solution of \3.1\) . 

A few remarks are now in order. 

Remark 3.7. It is clear that the eigenvectors in (|3.15p satisfy that e e — > eo, 
e* — > eo when e — > 0, where eo is as in (|3.17p . Therefore, by (|3. 18|) we have 
6(e) —> as e —> 0. 

Remark 3.8. For a given e > 0, if 6(e) 7^ 0, then the equation (|3.ip 
has a saddle-node bifurcation at A* < Aq as shown in Figure 13.1( a): if 
6(e) = 0, then (|3.ip has an (supercritical) attractor bifurcation as shown in 
Figure l3~lT b) . 
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Figure 3.1. (a) Saddle-node bifurcation and hysteresis, and 
(b) supercritical attractor bifurcation. 



To apply the above dynamic bifurcation theorem it is crucial to verify 
the global asymptotic stability for (|3.16p at A = Ao- The following theorem 
plays an important role to ensure this condition, which was proved in [9"llll|. 

Theorem 3.9. Under the conditions \3. 5\) , \3. 6\) and \3.18\) . ifG(e ,X ) / 
0, then u = is a globally asymptotically stable equilibrium point of li3.16\) 
at A = Aq. 



4. Dynamic bifurcation and stability of the Couette-Poiseuille 

Flow 

4.1. Functional setting. We recall here the functional setting and some 
basic mathematical properties of equations (|2,2ip with boundary conditions 
given by ()2.14|) . The spatial domain is M = (ri,ri + 1) x (0,L), the coor- 
dinate system is (r, z), the velocity field is u = (u z ,u r , ug) and u = (u z ,u r ). 
Let 



(4.1) 



H = {u= (u,u e ) G L 2 (M) 3 
V = {u£H 1 (M) 3 C\H\ u = 



dWu = 
at r 



Hi = {ue H 2 (M) 3 n V | ?t = ^ 



0, u-n\ dM = 0}, 
= ri,r 2 , u z = at z 
= at z = 0,L}, 



0,L}, 



where H k (M) is the usual Sobolev spaces. 

The equations (I2.2ip are two-dimensional, therefore the results concern- 
ing the existence and regularity of a solution for (|2.2ip with (|2.14p are clas- 
sical. For each initial value uq £ H, (|2.2ip and (I2.23P with (|2.14p has a weak 
solution 

u G L°°([0, T],H) H L 2 ([0, T],V), V T > 0. 
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u£ C([0,T],y)nL 2 ([0,T], J ff 1 



VT>0. 



For every uq € H and & > 1 there exists a tq > such that the solution 
u(t,u ) of (123TD . OT23D with (|2TT4l satisfies 



u(t,u ) S H k (My 



V t>T . 



Furthermore, for each k > 1 and no € H k (M) s , there exists a number (7 > 
depending on k and </? which bounds the solution u(t, uo) of (|2.2ip and (|2.23[) 
with (|2.14p in the i? fc -norm (see 0): 



\u(t,u )\\ H k < C, 



V f > 0. 



Finally, it is essentially known that for any A, 7 and u the equations (|2.21|) 
with (12.141) have a global attractor; see Temam [14] and Foias et al. [2]. 



4.2. Eigenvalue problem of the linearized equations. The linearized 
equations of <|2.21 j> with (|2.14|) are given by 

( QiXi dj) 

- Au z - 7(1 - 2r)u r + 7(^0 + r(l - r))-^- + -£ = 0, 



(4.2) 



Au r + A(l - fi)ru e + 7 (W + r(l - r))^ + ^ = Au e , 



A^- 7 (^ + r(l-r))^ 



Aii r 



<9u z <9u r 



0, 



with the boundary conditions 
u = 



at r = n, ri + 1, 



(4.3) 



du r dug 
u z = 0, -t— = - T — = at z = 0, L, 

oz oz 



where r = r — 77, A = x/T, T is the Taylor number, T and 7 are given by 
(PT22jh and u is defined by (pET) . We denote 

e = (7, 1 — /j) e (0, 00) x (0, 00). 



Let Aq > be the first eigenvalue of (14.20 . (|4.3j) . and we call 



(4.4) 

the critical Taylor number. 



e\2 



(A ) 
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When e — > 0, the system (|4.2[) reduces to the following symmetric linear 
equations 

- Au z + / = 0, 

oz 



(4.5) 



Au r + — = Aug, 
or 

Auq = Xu r , 



0. 



du z du r 
dz dr 

It is known that the eigenvalue problem (j4.5j) with boundary conditions 
4.3j) has the eigenvectors as follows 



(4.6) u Q = (u z ,u r ,u e ) 



11-1 

(— sin azR'Or), — sin azR f (r), cosaz 
a a a Aq 



dz 2 



a 2 ) i?(r)), 



where a is the wave length which is given by (see [T]) 

Kir 
a = — ' 

for some if > 1, with if depending on L. When L is large, if is taken to 
be the minimizer of 



mm 

k 



^-3.117 



The number Ao is the first eigenvalue of (|4.5|) . and when a ~ 3.117, Aq 
1700. 

Consider 



a 1 R 



dr 2 

R = 0, r! = o, 



-Q 2 Aq-R, 



(ir 2 



i? = at r = n, ri + 1. 



Its solution can be expressed as (see Chapter 11-15 of pQ) 
(4.7) R(r) = cos otox — Pi cosh u\x + cos + fa sinh aix sin 
where x = r 1 77 



(4.8) 



2 ,r = r-n, 

0i = 0.06151664, fo = 0.10388700, 
a = 3.973639, a x = 5.195214, a 2 = 2.126096. 

In addition, the first eigenvalue Ao of (|4.5p with (|4.3p is simple except for 
the following values of L: 

kir 



(4.9) 



U 



a 



k = 1,2,- 



a ^ 4.2. 
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4.3. Secondary flows and their stability. We are now in position to 
state and prove the main theorem in this section. We remark here that 
under conditions (12.111) and (12. 15ft with scaling d = 1, the condition /x — > 1 
can be replaced by 

(4.10) (1 - n) ri = 2 + 8 

for some 8 > and r\ sufficiently large. This condition implies that 

a = (r ] 2 - fi)/(l - rf) ~ 8. 

Theorem 4.1. Assume that \4- Mty holds true, and L ^ is sufficiently 
large, where is given by Then there are R > and 8± > such that 

for any r\ > R and \^\ < 8\, the following assertions hold true, provided 
that the Taylor number T = A 2 satisfies that < T — T c < 82 with T c as in 
f or some 82 > 0, or equivalently < A — Xq < 82 for some 82 > 0: 
(1) There exists an attractor T, £ x = {u^u^} C H of 12.21]) with {2.14% 
where u\ and are steady state solutions given by 

u\ = ai(X,e)u + vi(X,e), 
U2 = -cr 2 (A, e)u + v 2 (X, e), 
k Vi{X,e) = o(\ai(X,e)\), i = 1,2. 
Here uq = (u z ,u r ,ug) is the eigenfunction given by {4-oV , and 



(4.11) 



(4.12) 



^/V\e)+A(3 £ {X) + \b{e)\ 
*i,2 = ~ 2C + o(|6|, |/?!|), 

/3f(A) = a(e)(X - X E ) + o(\X - Ag | ) , a(0) > 0, 

b{e) = - [ {w £ ■ V)w £ ■ w £ *dx, 
Jm 

where w £ = {w £ ,w £ d ), with w £ = (w £ z ,w £ ), is the eigenfunction of 
and j^.3[ ) corresponding to Aq, w £ * is the dual eigenfunction, 
and C > a constant. 

The attractor S| attracts all bounded open sets of H \ T, where V 
is the stable manifold of u = with co-dimension one. Namely, the 
space H can be decomposed into two open sets and U^: 

H = u\ + u^, t/^n 1/2=0, e r = du? n du£, 

with G (i = 1,2), such that 

lim \\u(t, (p) — u^\\ L 2 = 0, when ip G (i = 1, 2), 

t^oo 

where u(t, (p) is the solution of 112.21]) with \2.13\) and §2. 



Remark 4.2. As mentioned in the Introduction, the existence of a bifurca- 
tion to a steady state solution is classical, and can be proved by the classical 
Krasnoselskii theorem. There are two new ingredients in this theorem. First, 
here it is proved that the the basin of attraction of the bifurcated attractor 
S| is H\T. Second the explicit form of the bifurcated solutions. These new 
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ingredients enable us to determine, in the next two sections, the asymptotic 
structure in the physical space of the solutions of the problem, leading to 
justification of the flow structure of the TCP flow. 

Proof of Theorem \4-l\ We shall apply Theorem 13.61 to prove this theorem. 
Let H and H\ be defined by (14, lh . We define the mappings 

f L A = A + XB : Hi -» H, 

L\ = L\ + C\ : Hi —> H, 



(4.13) 
by 



Au = {PAu, Au e }, 
Bu = {P(0,u e ),u r }, 

C e x u = |p^ 7 (l-2rK + 7(Wo + r(l-r))^, X{l-fj,)ru e 

+ 7(^0 + r(l - r(l - f))^\ , 7(^0 + f(l - 

G{u) = {-P(u-V)u,-(u-V)u e } . 

where u = (u z ,u r ,uo) G JETi, and the operator P : L 2 (M) 3 — > i/ is the 
Leray projection. Thus the equation (|2.2ip with (|2.15|) can be written in 
the abstract form 

= L x u + C\u + G{u). 

It is well known that the above defined operators satisfy the conditions 
([32])- ([331) • In particular 

B : H a — > H is bounded for all a > 0, 
C\ : H a — > H is bounded for all a > — , 

G : H a — > H is analytic for a > — , 

where H a is the fractional Sobolev spaces defined by the interpolation be- 
tween Hq = H and H\. 

It is clear that L\ is symmetric, and G : H a — > P (a > ^) is bilinear 
satisfying that 

(G(u,v),v)h = 
Hence, the condition (|3.18p holds true. 

With the operators in (|4.13j) . the eigenvalue equation (|3.17p corresponds 
to the problem (|4.5p with (|4.3p , and the eigenvector eo is precisely the uq = 
(u z ,u r ,ug) given by (|4.6p . 

Now consider the eigenvalue problem 

(4.14) L\u = /3(X)u, u = (u z ,u r ,u 9 ) € Hi. 
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By (I4.13p . the abstract form (j4.14j) is equivalent to the following eigenvalue 
equations in H\: 



(4.15) 



Au r 



dp 
Or 



+ Xug = f3(\)u r 



Aug + Xu r = (3(X)ug, 

divn = 0. 



It is known that the eigenvalues (5k (k = 1, 2, • • • ) of (|4. 15)) in H\ are real 
numbers satisfying 

'ft (A) > /?2(A) > ••• > PkW >•••, 
Pk — ► -co as k —>■ oo. 



(4.16) 



Because the first eigenvalue Ao of (|4.5p with (|4.3p is simple, the first eigen- 
value /3i(A) of (|4.15p with (|4.3p is also simple at A = Ao- Noting that for the 
first eigenfunction no of (|4.15p with (|4.3p at A = Ao, which is given by (|4.6j) . 
we have 



{Bu ,u ) H 



2 / u r ugdx 

JM 



a 2 A 



o Jo 



a 2 A 



cos 2 azdz 



cos 2 azdz 



n+l 



i? 



'I 



70 



n+l 



dz 2 
dz 2 



2 



a 2 ) i* 



fir 



> 0. 



Hence, from Theorem 13.31 and (|4.16p we can derive the conditions (|3.5p and 
(1331) at A = A . 

Finally, by direct computation, we have 

G(u ) ± 0, 

for the first eigenfunction no of (|4.5|) with (|4.3p . By Theorem 13.91 one can 
derive that u = is a globally asymptotically stable steady state solution of 
the following equations with (|4.3I) at A = Aq: 



dp 



Auz ~ if 
dz 



(u ■ V)u z 



du z 
dt 

du r dp _ 

— = Au r + Xug - — - (n • V)n r 

<9n< 
~dt 

divn = 



9r 

Aug + Au r — (n • V)ug, 



which correspond to the abstract equation (|3.16p . Thus, this theorem follows 
from Theorem 13.61 The proof is complete. 
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5. Structural Transition of the Couette-Poiseuille Flow: 
Boundary-Layer Separation 

In this and the next sections, we study the structure and its transitions in 
the CP and TCP flows, using a recently developed geometric theory of 2D 
incompressible flows by the authors; see a recent monograph by the authors 
[12j and the references therein. The results obtained provide a rigorous 
characterization on how, when and where the Taylor vortices are generated. 

5.1. Geometric Theory of Incompressible Flows. We first give a re- 
capitulation of some results in the geometric theory of incompressible flows, 
which will be used in characterizing the structure and its transitions in the 
CP and TCP flows. We refer interested readers to above references for 
further details of the theory. 

5.1.1. Structural stability. Let M C R 2 be an open set, and C r (M, R 2 ) be 
the space of all C r (r > 1) vector field on M. We denote 

D r (M,R 2 ) = {v G C r (M,R 2 ) | divv = 0, v ■ n\ 9n = 0} , 
B r (M, R 2 ) = Ive D r (M, R 2 ) | ^ = on dSl 

[ OT 

Bq(M, R 2 ) = {v G D r (M,R 2 ) | v = on dU} , 

where r, n are the unit tangent and normal vectors on dM. 

Let X = D r (M, R 2 ), or B r (M,R 2 ), or B^M^ 2 ) in the following defini- 
tions. 

Definition 5.1. Two vector fields u, v £ X are called topologically equivalent 
if there exists a homeomorphism of (p : M — > M , which takes the orbits of u 
to orbits of v and preserves their orientation. 

Definition 5.2. A vector field u £ X is called structurally stable in X if 
there exists a neighborhood O C X of u such that for any v G O, u and v 
are topologically equivalent. 

For u G Bq(M, R 2 ) (r > 2), a different singularity concept for points on 
the boundary was introduced in |7j, we proceed as follows. 

(1) A point p G dM is called a 9-regular point of u if 

du T (p) ^ Q 
dn ' 

and otherwise p G dM is called a 9-singular point. 

(2) A (^-singular point p G dM of u is called nondegenerate if 



det 



/ d 2 U T ( P ) d 2 U T (p) ' 

drdn dn 2 



^0. 

A nondegenerate 3-singular point of u is also called a 9-saddle point. 



i d 2 u n (p) d 2 u n (p) 
\ drdn dn' 2 
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(3) A vector field u £ Bq(M, M 2 ) (r > 2) is called -D-regular if it is regular 
in M, and all d-singular points of u on dM are nondegenerate. 
The following theorem provides necessary and sufficient conditions for 
structural stability in Bq(M, M 2 ). For the structural stability theorems in 
D r (M, M 2 ) and B r (M,R 2 ), see [121 [8]. 

Theorem 5.3 (Ma and Wang [7]). Let u £ £g(M,R 2 ) (r > 2). T/ien it is 
structurally stable in Bq(A!,M. 2 ) if and only if 

(1) u is D-regular, 

(2) all interior saddle points of u are self-connected; and 

(3) each d-saddle point of u on dM is connected to a d-saddle point on 
the same connected component of dM . 

Moreover, the set of all structurally stable vector fields is open and dense in 
Bq(M, M. 2 ). 

5.1.2. Structural bifurcation and boundary-layer separation. Let u(-,X) £ 
Bq(M,M. 2 ) be a one-parameter family of vector fields. Assume that the 
boundary dM contains a flat part T C dM and xq £ T. For simplicity, we 
take a coordinate system (xi,X2) with xo at the origin and with T given by 

T = {(xi,0) \\xi\ < 6 for some 5 > 0} . 

Obviously, the tangent and normal vectors on T are the unit vectors in the 
x\- and X2-directions respectively. In a neighborhood of xo £ T, u(x, A) can 
be expressed near x = by 

u(x, A) = X2v(x, A). 
To proceed, we consider the Taylor expansion of u(x, A) at A = Ao: 
U(x,A) = n°(x) + (A-A )n 1 (x)+ O (|A-A |), 
\ u°(x) = x 2 v°(x). 

Let u° = (115,1^2). We assume that 

(5.2) ^-AV-O. 

(5.3) ind(u°,0) = 0, 

(5.4) M5)^0„ 

(5.5) -r U ; / 0, for some k > 2. 

a K Tdn 

The following is a theorem describing the structural bifurcation on bound- 
ary for the case with index zero, which characterizes the boundary-layer 
separation for 2-D incompressible fluid flow. 

Theorem 5.4. Let u(-,X) £ Bq(M, M 2 ) be as given in 115.1]) satisfying the 
conditions \5.ty) - l5~5\) . Then the following assertions hold true: 
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(1) As X < Ao (or A > AoJ, the flow described by u(x, A) is topologically 
equivalent to a parallel flow near xq S V, as shown in Figure \5.1h . 

(2) As Ao < A (or X> Xq), there are some closed orbits of u separated 
from xo £ T, i.e. some vortices separated from xq £ V, as shown 
schematically in either Figure HPl (c) or (d). 

(3) If k = 2 in $5. 5|) . the vortex separated from xq S T is unique, and 
the flow structure enjoys the following properties: 

(a) there are exactly two d-saddle points x^ = (xf, 0) G T ofu(x, X) 
near xq = with x^ < < x^, which are connected by an 
interior orbit 7(A) of u(x, A); 

(b) the closed orbits separated from xq are enclosed by the interior 
orbit 7(A) and the portion of the boundary between x~ and x + ; 
and 

(c) the interior orbit 7(A) shrinks to xq as X — > Aq. 




(O (d) 

Figure 5.1. Boundary-layer separation and re-attachment 
near a flat boundary. 

We mention that when dM is curved at xq G dM, corresponding theorem 
is also true; see [121 El H] . 

5.2. Boundary-layer separation in the TCP flow. Hereafter, we al- 
ways assume that the conditions in Theorem 14.11 hold true. The main objec- 
tive of this and the next sections is to study the asymptotic structure and 
its transition in the physical space in the CP flow and the TCP flow regimes. 
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As mentioned in the Introduction, the results obtained in these two sections 
provide a rigorous characterization on how, when and where the Taylor vor- 
tices are generated. In particular, contrary to what is commonly believed, 
we show that the propagating Taylor vortices (PTV) do not appear after 
the first dynamical bifurcation, and they appear only after further increase 
the Taylor number so that the boundary-layer or interior separation occurs. 

We shall prove that the type of separations (boundary-layer or interior) is 
dictated by the (vertical) through flow driven by a pressure gradient along 
the rotating axis. Hence we consider two cases. The first is the case where 
through flow vanishes at the cylinder walls, i.e. Wo = in (|2.7p . This leads 
to boundary-layer separation, and will be studied in this section. 

The second is the case where Wo 7^ 0, leading to interior separations. 
This case will be addressed in the next section. 

5.2.1. Main results. Consider now the case where Wo = 0, and the main 
results are the following two theorems. The first gives a precise characteri- 
zation on how, when and where the propagating Taylor vortices (PTV) are 
generated for the bifurcated solutions, and the second theorem describes 
time evolution to the PTV for the solutions of the equations. 



Theorem 5.5. Assume that the conditions in Theorem \4-l\ hold true, and 
assume that the constant velocity in boundary condition {2. 7| ) vanishes, i.e., 
Wo = 0. Then there is a 70 > such that for any < 7 < 70 where 7 is 
defined by \2. 22\) , the following assertions hold true for the two (bifurcated) 
steady state solutions v} = U cp + u} (i = 1, 2) of (2lfy and pTTp , where U cp 
and are given by \2.3\) and ftJTTTty: 

(1) For v^, there is a number Ai > Aq (Aq as in ^4-4] )) such that 

(a) for Aq < A < Ai, the vector field Vi = (W + u\ z ,u\ r ) is topo- 
logically equivalent to the vertical shear flow (W, 0) as shown in 
Figure WM, and 

(b) for Ai < A, there is a unique center of v^ separated from each 
point (z n ,ri) £ dM, where z n ~ (4ra + l)-7r/2a (resp. {z n ,r<2) 6 
dM, z n = (4n + 3)ir/2a) for n = 0, 1, • • • ,ko, as shown in 
Figure [5731 with a as in J^.ffi ). 

(2) For V2, there is a number bo > such that only one of the following 
two assertions holds true: 

(a) // \ b(e)\ < bo (b(e) as in \4.12ty ), there is a A2 > A§ such that 



(i) for Aq < A < A2, the vector field = (W + u^ z , u^) has 
the topological structure near r = r\ (resp. r = r2) as 
shown in Figure Ejj and 

(ii) for A2 < A, there is a unique center ofv\ separated from 
each point (z n ,ri) £ dM (resp. (z n ,r2) G dM) for n = 
0, 1, ■ • • , ko, as shown in Figure [5751 

(b) // \b(e)\ > bo, then for A > A2 = Ag ; the vector field is 
topologically equivalent to the structure as shown in Figure \B~4\ 
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(3) There exists a A3 > Ai (Ai > A2) such that for Aj < A < A3, the 
vector field (i = 1,2) is topologically equivalent to the structure 
as shown in Figure and A3 is independent 0/7. 



From Theorems 14.11 15.51 and 15.31 we immediately derive the following 
theorem, which links the dynamics with the structure in the physical space. 



Theorem 5.6. Assume that the conditions in Theorem \4-l\ hold true, Wq = 
in {2.1), and < 7 < 70. Then there are Xj (1 < j < 3) with Ai = Ai, 



A* < A2 < A2, A3 = A3, where Xj (1 < j < 3) are as in Theorem 15.51 and 
X* is the saddle-node bifurcation point (if it exists), otherwise X* = A§, and 
for each A < A3 the space H can be decomposed into two open sets 

H = ul + ui, n C/ 2 A = 0, 

such that the following assertions hold true. 

(1) For ip G J7- there is a time to > such that if X < Aj (i = 1,2) 
/or i/ie solution u(t,tp) = (u,u s ) of i2.2\) with fg.Tp and fff.&g)) , t/ie 
vector field u is topologically equivalent to the structure as shown in 
Figure [57B. and if Aj < A < A3, u is topologically equivalent to the 
structure as shown in Figure \5.4\ for all t > to • 

(2) In particular, for the initial value ip G near the Couette-Poiseuille 

flow \2. 3\) . if X > Aj there is a time to > such that the solution 
u(t, if) of 112.2]) with {2.7\j and \2.23\) has a boundary-layer separation 
at t = to, i.e., the vector field u ofu(t,(p) is topologically equivalent 
to the structure as shown in Figure HOI for t < to, and u is topo- 
logically equivalent to the structure near boundary dM as shown in 
Figure [Ql for t <t. 



r 

r=r+l 



z 



r=r i 



Figure 5.2. Vertical shear flow with vanishing boundary velocity 
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Figure 5.3. Boundary separation and re-attachment. 




Figure 5.4. Propagating Taylor Vortices after boundary- 
layer separations. 



Remark 5.7. The critical values Xi (i = 1,2) in Theorem 15.61 which deter- 
mine the boundary-layer separation for the Taylor-Couette-Poiseuille flow 
problem depend continuously on the non-dimensional parameter 7 = p$d? / 'Apu 2 . 
The expansion is as follows 

(5.6) ^1(7) = a il + °(ItI ) ? ai > a constant, i = 1, 2. 

In addition, the condition < 7 < 70 in the above theorems is the conse- 
quence of the narrow-gap assumption (i.e. d is small). 

Remark 5.8. The height L has an integer number times the length of it /a, 
i.e., L = K ■ — for some integer K. Hence the number of vortices separated 
= 7*1 and r = r2 in Theorems [ 

K + 1" 



from r 



kn 



I and 15.61 satisfies 
for r = n (or r 



r '2) 



K 
~2 



for r = r2 (or r = r%) 
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where [a] is the integer part of n. 

5.2.2. Proof of Theorem 15. 51 We proceed in a few steps as follows. 

Step 1. By Theorem l4.lt as the Taylor number T = A 2 satisfies T c < T 
(T c = Xf), or equivalent^ X e < A, the equations (pT2j) with (pTTj) and (pT4l 
generate from the basic flow (|2.3|) two steady state solutions 



(w(r) + ui,ui,V(r)+u^j, * = 1,2, 



which are asymptotically stable under the axisymmetric perturbation. 
It is clear that the number b(e) in ()4.12|) satisfies 

Iim6(e) = 0, (e = (l-n,j)). 

Hence, to prove this theorem, it suffices to proceed only for the following 
two vector fields: 



(5.7, vt = iW(r) , 0) + ^±^EM^^ 



(5.8) ^. (H r (r) ,o ) + ^±^M±» Sci 

where b = b(e), a = a{e) > near e = 0, C > a constant, 
W(r) = —jr(l — r), (r = r — ri), 

So = f — sin azR'(r), — cos a-Z-R(r) ] , 
V a / 

and R(r) is given by (|4.7p . 

Step 2. Claim: V± has a boundary-layer separation on Ai < A /or some 
Ai>A§. 

We shall apply the structural bifurcation theorem, Theorem 15.41 to prove 
this claim. To this end, let 



A = (Jb 2 + 4a(A - A§) - 6^ /2aC, 



which is an increasing function of A with A = at A = Aq. Thus the vector 
field (15.71) can be rewritten as 



(5.9) Vi = (A sin azR'(r) - jf(l - r), -aAcos azR(r)) . 



We discuss the structural bifurcation of (|5.9p only on the portion r = r\ 
of OM, and on this portion r = r*i + 1 of dM, the proof is the same. 

Since R(r\) = R'(ri) = 0, the zero points of = on r = r\ are 
given by the solutions of the equation 

(5.10) A sin azRl'in) - 7 = 0. 



(5.12) 
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Let Ao = r y/R"(ri). Then, we infer from ()5.10p that for r = r%, 

d yx { ^0 whenA<A , 

(5-11) -ttM n 1 . . (4n + l)7r 

dn = when A = A and z n = - — . 

v 2a 

As in dEU, 

can be Taylor expanded at A = Ao as follows: 

'y 1 A = yo + (A-A )F 1 , 
y° = y r °) = (A sin azR' - 7f (1 - f), -aA cos azR), 
V 1 = (V},V r ) = (smazR', — acos azR). 
We derive from (|4.7p that 

(5.13) R"(r) = — ag cos ao^ 

+ [/3i(«2 ~~ a i) + 2/?2«ia2] cosh(aix) cos(a2^) 
+ [foi&i — ot\) + 2/3iai«2] smh(a<ix) sin(«2a;), 
where x = r — ^ . Inserting (|4.8p into f)5 . 13|) we obtain 

(5.14) R"(n) ~ 8 + 1.7ei ~ 28. 
Thus we derive from (|5TT]) . (f5T2l) and (I5T41 that 

<9F° <9V° , . /4n+l 



9n <9r 2a 
dV 1 



dn dr 



o at (z n ,n) = ( n ^ 7r,n 

(28,0) at(z n ,n), 
28a 2 Aq at (z n , ri). 



dr 2 dn dz 2 dr 

Namely, the conditions (|5.2p . (|5.4p and (|5.5p are valid. 

Now, we verify the condition (|5.3p . i.e., we need to prove that 



(5.15) ind^^,(^,ri)J = 0. 

To this end, we consider the equation 
dV x 

(5.16) — — = A sin azR"(r) - j(l - 2r) = 0. 

ar 

The function R"(r) has the Taylor expansion at r = r\ as follows 

(5.17) R"{r) = R"(n) + i?"(n)r + o(\f\). 
From (|4.7p and (|4.3p we can derive that 

17 5 

(5.18) ft w (ri) ~ -64 x — - 10 x el ~ -175. 
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On the other hand, we note that 

7A 7A 



A 



R"(n)A 28A 



Thus, by (|5.14p . (|5.17p and (|5.18p . the equation (|5,16p near (z,r) = (z n ,ri] 
can be rewritten in the following form 



/A \ / 175Asinaz \ „ ,. 

(5.19) [l - - sin a, J + {-^~ ~ *) r + o(\f\) = 0. 



Obviously, there is a 5 > such that the equation (|5.19p has no solutions in 
0< r — n = r < 5, \z — z n \ < 5 and A < An. Namely, in a neighborhood of 

(z n ,ri), the vector field has no singular points for A < Ao. Hence, by 
the invariance of the index sum of singular points in a domain, one can get 
(J5TT5]) . 

Let Ai be the number that A(Ai) = Ao. It is clear that Ai > Ag. Thus, 
this claim follows from Theorem [57 



Step 2. The vector field has no singular points in M for Aq < A < Ai, 
and has only one singular point in each domain (2^,^4-1) x (ri,r± + 1), 
z k = !f (0 < k < K = - 1), for all A > Ai. 

From (|5.9p we find that the singular points of must be in the lines 
4 = {(^ + £,r) I n < r < n + 1} C M, and 



^ir = 0, ^ = (-1)^(0 - 7 r(l-r), on 4 

Therefore, to prove this claim we only consider this equation: 

, t R'(r) An . 
(5-20) (_i)*_^__o ? .(i- ? .) = o, 0<r<l. 



It is known that i?'(r) is odd and r(l — r) is an even function on the variable 
x = r — I . These functions are illustrated by Figure 15.51 



SEPARATIONS IN TAYLOR-COUETTE-POISEUILLE FLOW 



27 




r (1-r ) 



1 



r 



Figure 5.5. r is a zero point of (|5.20p . 



Direct computations show that (|5.20[) has no zero points in < r < 1 
for A < Ao, and there is a unique zero point Tq G (0,1) for all Ao < A. 
Furthermore, < Tq < \ as k = even, and ^ < Tq < 1 as k = odd. Thus, 
this claim is proved. 

Step 3. By the invariance of index sum of singular points in a domain, it 
is easy to see that the singular point of in the domain (z^, Zk+x) x (ri, ^2) 
must be a center, which is enclosed by an orbit of connected to both 
boundary saddle points Z\ and Z2 , as shown in Figure 15.31 ?L > ?2 are on 
r = n as k = even, and on r = r<i as k = odd. 

Finally, by Steps 1-3 one readily derives this theorem. The proof is com- 
plete. 

6. Structural Transition of the Couette-Poiseuille Flow: 



In this section, we study the structure and its transitions in the TCP 
flow for the case where the z-direction boundary velocity Wo ^ in f|2.7j) . 
We show that the secondary flow from the Couette-Poiseuille flow (|2.3p 
will have interior separations as the Taylor number T exceeds some critical 
value T which is an increasing function of the nondimensional parameter 
7 = pod s /4pi>. For this purpose, we recall some rigorous analysis by the 
authors on interior separations of incompressible flows; see |10t 112] . 

6.1. Interior structural bifurcation of 2D incompressible Flows. 

Following notations used in Section [5.11 we let u(-,\) G D r (M, M 2 ) (r > 1) 
have the Taylor expansion at A = Aq as follows 



Interior Separations 



(6.1) 



u[x, A) = u\x) + (A - AoK(x) + o(|A - A |). 
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We assume that xq £ M is an isolated singular point of uq(x), and 



(6.2) 
(6.3) 
(6.4) 



ind(ii°, xq) = 
Du°(x ) 
v^ixo) • e 2 / 0, 



f dul(x ) du°(xo) \ 

dxi 8x2 

i du^(x ) dv% {x ) - 

\ dxi dx2 / 



where e 2 is a unit vector satisfying 



(6.5) 



Du°(xo)e2 = ae\ 
Du°(x )e 1 = 0. 



(«/0), 



We also assume that u° £ C m at xq £ M for some even number m > 2, and 



(6.6) 



d k (u°(x ) • e 2 ) 



0, 1 < k < m = even, 
7^ 0, k = m = even. 



Then we have the following interior structural bifurcation theorem, which 
was proved in [lQl 02] . 



Theorem 6.1. Let u(-,A) £ D r (M, R 2 ) be as given in $6.1\) satisfying the 
conditions Ifi-j^ - pT^} ), T/ien i/ie following assertions hold true. 

(1) As A < Ao (or A > Aoj, i/ie /Zow described by u(x, A) is topologically 
equivalent to a tubular flow near xq £ M as shown in Fiaure [b\lV a). 

(2) As A > Ao (or A < Xq) there must be some centers ofu(x,X) sepa- 
rated from xq £ M as shown schematically in either Figure ROl (c) 
or (d). 

(3) The centers are enclosed by an extended orbit 7(A), and 7(A) shrinks 
to xq as A — » Ao- 

(4) // i/ie condition \6. 6\) is satisfied, then the center separated from 
xq £ M is unique, as shown in Figure UL~T\ (c). 
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(c) (C) 

Figure 6.1. Interior separation. 

6.2. Interior separation of the Taylor-Couette-Poiseuille Flow. 

6.2.1. Main theorems. Let the constant velocity W$ ^ in the boundary 
condition (|2,7p . Then, we have the following interior separation theorems 
for the TCP flow problem. The results obtained here are in agreement with 
the numerical results obtained by Raguin and Georgiadis; see e.g. Figure 8 
in p3]. 

Theorem 6.2. Assume that the conditions in Theorem \4-l\ hold true, and 
let Wo ^ in {2. 7| ). Then there exists 70 > such that if < 7 < 70, then 
for the two (bifurcated) steady state solutions = lA cp + (i = 1,2) of 
\2.2\) and \2. 7| j, the following assertions hold true: 

(1) For Vi, there is a \\ > \ e such that 

(a) for Aq < A < Ai, the vector field v\ = (W + u\ z , u\ r ) is topolog- 
ically equivalent to the structure as shown in Figure Xb\2\ and 

(b) for Ai < A, there is exactly a pair of center and saddle points 
separated from a point in each domain (4, 4+i) x ( r i> r 2) C 
M, with Zk = ^ (0 < k < K), as shown in Figure [thW c). 
Moreover, the value X\ is a continuous and increasing function 
of J- 
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(2) For t>2 there is a number bo > such that only one of the following 
two assertions holds true. 

(a) If 1 6(e) | < bo, where 6(e) is given by \4.1%fy ), then there is a 
A2 > Aq (^i — ^2) such that the same conclusions as Assertion 
(1) holds true for v% with A2 replacing Ai . 

(b) If\ 6(e) I > 60, then when A > A2 = Aq, the vector field v% = (W+ 
u 2zi u 2r) i- s topologically equivalent to the structure as shown in 
Figure 

(3) There exists a A3 > Ai, with A3 independent 0/7, such that if A, < 
A < A3, the vector field (i = 1,2) is topologically equivalent to the 
structure as shown in Figure [g73i 



The following theorem is a direct corollary of Theorems 15.31 14.11 and 16.21 
and provides a link between the dynamics and the structure in the physical 
space. 



Theorem 6.3. Assume that the conditions in Theorem \4-l\ hold true, Wo 7^ 
in (2.7), and < 7 < 70. Then there are Aj (1 < j < 3) with X\ = X\, 



A* < A2 < A2, A3 = A3, and for each A < A3 the space H can be decomposed 
into two open sets 



H = u\ + U2, n C/ 2 A 



such that the following assertions hold true. 

(1) For ip G U£ there is a time to > such that when t > to, for the 
solution u(t, <p) = (u, ug) of H2.2\) with \2.7\) and H2.23\) . as A < Aj the 
vector field u is topologically equivalent to the structure as shown in 
Fiaure [672\ and as Aj < A < A3 (i = 1,2) u is topologically equivalent 
to the structure as shown in Figure [TOl 

(2) In particular, for the initial value cp £ near the Couette-Poiseuille 
flow \2.ty) . and for A > Aj (i = 1, 2), there is a time to > such that 
for the solution u(t,<p) = (u,u$) of \2. 2\) with {2.7\) and k2.23\) . u 



has an interior separation at t = to from a point in each domain 
(5fc,Zfc+i) x (^1,^2) C M, where Zk = ^f, as described by Theo- 
rem 1 6. 2\ with t replacing A. 
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Figure 6.2. Vertical shear flow with a constant velocity 
Wq / on boundary. 




Figure 6.3. Propagating Taylor vortices separated from the interior. 



6.2.2. Proof of Theorem 1 6. This proof is similar to that of Theorem 15.5 
Here, we only need to prove the interior structural bifurcation in Asser- 
tion (1) for the following vector field 

(6.7) V? = (A sin azR'(r) - jr(l -r) — jWq, —aA cos azR(r)) . 

We proceed by applying Theorem 16. 11 Consider the equations 
f ±A#'(r)- 7 r(l-r)- 7 Wo = 0, 
(6 ' 8) [ ± AR"(r) - 7 (1 - 2f) = 0. 



By computing, we can know that the equations (|6.8p have a unique solution 
(A Q ,rf ) with r£-n 



the maximum points of R'(r), i.e. 
(6.9) R"(rt) = 0, 



1 — (r — n) and < r+ < < r Q , where r J are 



l „- 
5' 
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The vector field (|6.7f) has the following expression at A = Ao 

(6.10) _ 
y 1 A = yo + (A-A )F 1 , 

V° = (F z °, V r °) = (A sinazi?' - 7 r(l - r) - jW , -aA cos azR), 
V 1 = (V*,Vr) = (sin azR', -a cos azR). 

It is easy to see that for the solution (Ao,r ) of (|6.8p the points (z2fc,?"o") 
and (^2fc+i) r cT) ^ M are singular points of V°, where z m = ^ + 

For simplicity, we only consider the structural bifurcation of V° at the 
singular point (z ,r ) = (^,rj). 

We can check that when A < Ao, the vector field V± given by (|6.7p has 
no singular point in M, and (£o> r o) is an isolated singular point of at 
A = Aq. Therefore we have 

(6.11) md(V°,(z ,r )) = 0. 
In addition, we see that 



(6.12) DV°(z ,r ) 



(dyl 9vl\ 

dz dr 

. dyl_ dyl_ . 

\ dz dr I 




(«0,n)) 

and the eigenvectors e\ and ei satisfying 

( DV (z ,r )e 2 = ae 1 , {a = a 2 AR(r )) , 
[DV o (z o ,r )e l = 0, 

are given by 

ei = (0,l), e 2 = (l,0). 

Hence, we have 

(6.13) V 1 (z ,r ) ■ e 2 = V z \z ,r ) = R'(r ) + 0. 

From (|6.11|) - (|6.13|) we find that the conditions (|6.2|) - (|6.4|) are satisfied by 
the vector field at A = A and (z,r) = (zQ,r ). 



Finally, we verify the condition (|6.6[) . By f)6.8|) we have 

d k V z °(z ,r ) 



.14) dk{V °-^ 



(*o,n>) d rk 



and 

(6.15) ^o,r o) =Ao ^ /(ro) + 27 



for k = 0, 1, 



By (|6.8p we see that 
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On the other hand, we can check that 

Ft"'(r) 2 

( 6 - 17 ) -WT>T-^ V Vn<r<r+ 

K"{r) 1 — 2(r — r\) 

where r+ satisfies (|6.9|) . From (|6.15|) - (|6.17|) it follows that 
( 6 - 18 ) + °" 

Thus, from (16.141) and (16.18h we derive the condition (16.61) with m = 2. 
By Theorem 16. 1\ the vector field has an interior separation from each 
point (y m ,r ) at A = A , where r = if y m = + ^ and r = if 
2/m = 71 ~ia' P ro °f is complete. □ 



References 

[1] S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, Dover Publica- 
tions, Inc., 1981. 

[2] C. FoiAS, O. Manley, and R. Temam, Attractors for the Benard problem: existence 
and physical bounds on their fractal dimension, Nonlinear Anal., 11 (1987), pp. 939- 
967. 

[3] M. Ghil, T. Ma, and S. Wang, Structural bifurcation of 2-D incompressible flows, 
Indiana Univ. Math. J., 50 (2001), pp. 159-180. Dedicated to Professors Ciprian Foias 
and Roger Temam (Bloomington, IN, 2000). 

[4] , Structural bifurcation of 2-D incompressible flows with the Dirichlet boundary 

conditions and applications to boundary layer separations, SIAM J. Applied Math., 
65:5 (2005), pp. 1576-1596. 

[5] S. Goldstein, The stability of viscous fluid flow between rotating cylinders, Proc. 
Camb. Phil. Soc, 33 (1937), pp. 41-61. 

[6] T. Kato, Perturbation theory for linear operators, Classics in Mathematics, Springer- 
Verlag, Berlin, 1995. Reprint of the 1980 edition. 

[7] T. Ma AND S. Wang, Structure of 2D incompressible flows with the Dirichlet bound- 
ary conditions, Discrete Contin. Dyn. Syst. Ser. B, 1 (2001), pp. 29-41. 

[8] , Structural classification and stability of incompressible vector fields, Physica 

D, 171 (2002), pp. 107-126. 

[9] , Dynamic bifurcation and stability in the Rayleigh-Benard convection, Commu- 
nication of Mathematical Sciences, 2:2 (2004), pp. 159-183. 

[10] , Interior structural bifurcation and separation of 2D incompressible flows, J. 

Math. Phys., 45 (2004), pp. 1762-1776. 

[11] , Bifurcation Theory and Applications, vol. 53 of World Scientific Series on 

Nonlinear Science, Series A, World Scientific, 2005. 
[12] , Geometric Theory of Incompressible Flows with Applications to Fluid Dy- 
namics, vol. 119 of Mathematical Surveys and Monographs, American Mathematical 
Society, Providence, RI, 2005. 
[13] L. G. Raguin and J. G. Georgiadis, Kinematics of the stationary helical vortex 

mode in taylorcouettepoiseuille flow, J. Fluid Mech., 516 (2004), p. 125154. 
[14] R. Temam, Infinite-dimensional dynamical systems in mechanics and physics, vol. 68 
of Applied Mathematical Sciences, Springer- Verlag, New York, second ed., 1997. 



34 MA AND WANG 

(TM) Department of Mathematics, Sichuan University, Chengdu, P. R. China 

(SW) Department of Mathematics, Indiana University, Bloomington, IN 
47405 

E-mail address: showEaig@indiaiia.edu 



